High momentum entanglement of cold atoms generated by a single photon scattering 
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With the mechanism of pairwise scattering of photons between two atoms, we propose a novel 
scheme to highly entangle the motional state between two ultracold neutral atoms by a single 
photon scattering and detection. Under certain conditions, it is shown that an arbitrary amount of 
entanglement can be obtained with this scheme. 
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Introduction. — Quantum entanglement, as one of the 
central topics in quantum mechanics, has been exten- 
sively studied on the Hilbert spaces with continuous vari- 
ables (CV) in recent years 0, i, i, i, |, 0, 0, Q. Com- 
pared to the finite dimensional entanglement on pho- 
tonic polarizations or atomic internal states, the CV sys- 
temprovides entanglement in a rich diversity of forms 
H B H H Q, which provide unique roles in quan- 
tum information processing [l|, [7j and in fundamental 
test for quantum physics [2]. It is only in CV systems 
that the unbounded high degree of entanglement be- 
comes possible, which has been studied recently for op- 
tical squeezed state optical parametric down conver- 
sion [3, 4], atom-photon scattering [5], and atom-atom 
motional state [6|, etc.. 

For motional entanglement between neutral atoms 
without interactions, it is verified that [6] one can not 
produce entanglement beyond 1 ebit with single-photon 
emission and detection, therefore, it unavoidably requires 
nonlinear bi-photon or multi-modes detections for high 
entanglement [6]. By introducing effective interactions 
between atoms in this letter, however, we find that it 
is possible to produce arbitrarily high entanglement by 
only a single photon scattering and detection, which is 
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FIG. 1: (Color online) Schematic plot of the physical model. 
The coupling photon is injected in a single-mode cavity along 
the y-axis, whereas the scattering photon is collected and de- 
tected along x-axis. The recoil of atoms is restricted in the 
x-direction. The inset shows the atomic configuration. 



more feasible and efficient for realistic implementations 
[7[. The basic physics for this scheme is that, for a cou- 
ple of cold atoms with the same electric dipoles, photon 
scattered from one atom can be efficiently re-absorbed 
and re-scattered by the other, this mechanism of "pair- 
wise scattering; of photons entangles the atomic 
wavepackets under the law of momentum conservation. 
Evaluated by the Schmidt number K 0, [H[ , it is shown 
that an arbitrary amount of entanglement can be pro- 
duced for ultracold atoms with proper physical control 
parameters. 

Theoretical analysis. — To describe the physical pro- 
cess, as shown in Fig. 1, it is assumed that a pair of 
identical cold atoms "a" and "b" are coupled to a sin- 
gle photon in a cavity along y-axis, with nearly paral- 
lel atomic electric dipoles d a>D , i-e., \d a \ = |d&| = d and 
their relative angle cp « 0. We assume the coupling pho- 
ton is scattered and therefore recoils the atoms to the 
x-direction where the photon detector is placed. Con- 
sidering the momentum exchange along x-direction, the 
Hamiltonian can be written as H total = + Hi , where 
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(1) 



is the x-dimensional momentum operator for the 
two-level atom "a" ( "b" ) with mass m and transition fre- 
quency LJ21, and o~ij = a) (a) is the creation (anni- 
hilation) operator for the coupling single mode with fre- 
quency uj c and coupling strength g c ; whereas a\ s (a^ s ) 
is for the scattering modes along x-axis with frequency 
ujk = c\k\ and polarization s and the coupling strength is 

9a(b)(k,s) ~- 



c\k\_ 



2eohV ' da(b) ' e k,s- 

We utilize the Dirac ket to denote the physical state, 
e.g., |l a , q a \ lb, qb\ k) represents that the atom "a" ("b") 
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FIG. 2: (Color online) Spatial probability distribution of the 
steady bipartite wavefunction D(x, k = k c ), with k c = 1, <r = 
0.2, (5 = 0.1 and E m — ^ [III]. The small figure shows the initial 
Gaussian distribution. 
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FIG. 3: (Color online) Density plot of momentum distribu- 
tions for the steady bipartite wavefunctions £>(q, k — kc) with 
k c = 1, £ = 0.1, and Em = (a) Different orders of 

scattered wavepackets are well separated with a = 0.3. (b) 
Wavepackets begin to overlap and interfere with a — 



has the internal state |l) a (6) and momentum hq a ^ : and 
the photon is scattered with momentum hk. In the 
slowly- varying frame, therefore, the system state can be 
expanded as: 



dq e 



-i(E/h+uj c )t 



-4(q)|ia,ga;26,g6;0) 



+B(q)| ;ih,%;0> + C(q)|l ,ga; 16,96; 0) 

+ ydqdfte- < ( B /' ,+w *> t £>(q ) A!)|l 0) g ;l 6) (a;fc), (2) 

where the bold symbol q simplifies the pairwise variables 
(q a ,qb), and the summation over the polarization s is 
implied in the integration over k. The kinetic energy is 
E = tf(ql + ql)/2m. 

In Eq. ©, D(q, fc) is the momentum wavefunction 
for the atoms and the scattered photon, which is our 
main concern for its induced atom-atom entanglement 
in the steady state. A(q) [-B(q)] is the instantaneous 
wavefunction for the atom "b" ("a") is excited, and we 
introduce M(q) = A(q) + B(q) for the symmetry. From 



the Schrodinger equation, the dynamical equations read: 
dM(q) 



dt 



dC(q) 
dt 

d£>(q, k) 
dt 



= -tAM(q) - 2ig* c C(q) - 7 M(q) 

- e iE ^' h Jdkj ds [g a (k)g b (k) x 
e i(^-c\k\)(t-s) M(qa+k ^ qb _ Ks) ^ 

= -ig c M(q), 

= — iexp [— i(u) c — c\k\)t] x 



(3) 
(4) 

(5) 

[g b (k)A(q a ,q b + k) + g a (k)B(q a + fc, ^)] , 



where the detuning A = — and 7 = 

27r Jcl/c |<7a(&)(fc)| 2 £(o;fc — ^c) is the atomic natural 
linewidth, and T = 2|g c | 2 7/A 2 denotes the scattering 
rate. E m = h 2 k 2 /m represents the kinetic energy mis- 
match generated when the coupling photon is scattered 
from one atom to another and be re-scattered back to the 
coupling mode @, [To[ . When this energy deficit is signifi- 
cant, the pairwise scattering happens only in a time scale 
r < h/E m and further cascaded processes will be elimi- 
nated [9]. In this work, we assume E m is small enough 
to allow needed cascaded pairwise scattering, and treat 
it as a constant 

The fourth term at the r.h.s. of Eq. ([3]) represents 
the pairwise scattering between two atoms. If the atomic 
dipoles are perpendicular (</? = 7r/2), this term vanishes 
and the pairwise scattering will be forbidden in the Eqs. 
([3]-[5]), which therefore reduce to the model of a single 
atom scattering [5]. When the dipoles are parallel ((^~0), 
however, photon scattered by one atom can be efficiently 
re-scattered by the other, therefore produces momentum 
correlation between them. 

Eqs. ([3]-[5]) can be analytically decoupled in the atomic 
position coordinates. We use the argument x= (x ai x b ) 
for wavefunctions to indicate their fourier counterparts, 
e.g., M(x) /dq e~* x>q M(q), where x • q = x a q a + x b q b . 
With weakly coupling conditions: g c ,j <C A, and 
E m /h <C T, the adiabatic solution of Eq. ([5j) yields: 
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[l_ e -^n(x,M)]/n( x ,^,t), (6) 
II(x,M) = i(u c -2\g c \ 2 //\-c\k\) 

+ r[l+ s'mc(E m t/h) cos(</?) cos (k c x a — k c x b )], 

where G(x) is the fourier counterpart of the ini- 
tial atomic momentum wavefunction G(q) = 
V^exp [— (q 2 + q b )/a 2 ]/ ^na which is generally specified 
as an unentangled Gaussian with momentum width 
a. For atoms with nearly parallel dipoles (cp ~ 0), 
we introduce cos(y?) = 1 — S 2 and use S (<C 1) as a 
dipole parallelity parameter in the following. N is the 
normalized factor and sinc(x) = sm(x)/x. 
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Generation of atom-atom entanglement with single- 
photon detection — Eq. (j6]) is the time evolution of collec- 
tive atoms-photon wavefunction, which will be projected 
to an atom-atom entangled state |6[ after the detection 
of the scattered photon. In the following, we neglect the 
light shift (2\g c \ 2 / A) and assume the photon clicks on a 
narrow band detector with wavevector k — +k c . 

When E m ^ 0, from Eq. ([6]), the steady atomic mo- 
mentum wavefunction reduces to a Bell-like state: 

L>(q, /c = /c c ,t^oc) oc G(q a + k c ,q b ) + G(q a , q b + k c ), (7) 

since for t ^> h/E m the energy-conserved single-atom- 
scattering process prevails |9| , which can produce at most 
1 ebit by a single-photon detection [6[. 

For E m = [ll|], however, the atom-atom pairwise 
scattering introduces effective interactions between them, 
which produces correlated wavepacket as: 

D(q, k = k c ,t^oo)(xD / (q a ^k c ,q b )^D / (q a ,q b ^k c ), (8) 
D (q) = exp [ — ] exp ( — ) x (9) 

°° 22 

2 exp [ ~^r (2n + 1)2] cos [^-( 2n+1 )(9»-*)i- 

n—0 c c 

Eq. ([9]) characterizes the essences for the pairwise scat- 
tering: the first Gaussian factor on the r.h.s. represents 
the bipartite total momentum which is provided by the 
initial momentum width cr; the Fourier series on the sec- 
ond line, which gives a periodic structure of the rela- 
tive momentum, characterizes different orders of pairwise 
scattering, with the total number of scattering orders 
controlled by the dipole parallelity parameter S through 
the second Gaussian factor on the first line. 

The spatial and momentum distributions of the atomic 
wavefunction of Eq. (|HJ) are plotted in Fig. 2 and Fig. 
3, respectively. In Fig. 2, one sees that, the scatter- 
ing of photon with wavevector k = k c correlates the 
atomic wavepackets with resonant spatial period A c , and 
produces grating-like wavefunction with highly localized 
atomic relative positions. Meanwhile, as in Fig. 3, the 
atomic momentum is correlated along q a +qb + k c = due 
to the momentum conservation. When the initial atomic 
momentum width a is well below the recoil momentum 
(a <C k c ) as in Fig. 3 (a), different orders of scattered 
wavepackets are clearly separated, which corresponds to 
a well entangled bipartite state; when a tends towards 
fc c , however, as in Fig. 3 (b), recoiled wavepackets begin 
to overlap and exhibit destructive quantum interferences 
for the high-order scattering, which will eventually de- 
stroy the momentum entanglement that can be seen more 
clearly in the following. 

The entanglement encoded in the wavefunction Eq. (|8|) 
can be quantitatively evaluated by the Schmidt decom- 



position of D(^k=k c ): 

CO 

D( q , k=k c ) = vXA(q), (10) 

n—— 00 

where S n (q) is the momentum Schmidt basis which is 
orthonormal and separable, with A n = 1. The degree 
of entanglement is then defined as the Schmidt number 

For the single-atom-scattering state in Eq. jTJ), com- 
pared with Eq. ([TO]) , one has at most K = 2 which is 
only 1 ebit [6[. For the pairwise-scattering-state in Eq. 
dH}, however, from the numerical results of K in Fig. 4 
in dependence of S and a, one sees that it may highly 
exceed the 1 ebit limit [6] for cold atoms (a < k c ) with 
parallel dipoles (J~0). 

For atoms cooled well below the recoil temperature, 
i.e., <j<C/c c , the Fourier series in Eq. (|9|) makes the rel- 
ative momentum highly localized in each single period, 
therefore, the summation over spatial frequency can be 
well approximated by a summation over discrete localized 
momentum modes. Along with Eqs. (|8 UTQ|) . it eventu- 
ally yields the Schmidt decomposition as a summation of 
different orders of scattered wavepackets: 

S n (q)=(-l) lnl+1 Sgn(n)G(q a -nk c ,q b +nk c +k c ),(ll) 
A n =7V'( e -^l"l-e-^l" +1 l) 2 , (12) 

where Sgn(n) is the signum function with value +1 for 
n > and — 1 for n < 0, and N f is a normalization factor. 
With Eq. ([T2]) . the degree of entanglement is obtained: 

K = l/Y,^l = yfi/S, (S«l). (13) 

n 

Eq. ([T3]) fits well with the numerical results when 
<j < k c /2 as shown in Fig. 4. It indicates that, by uti- 
lizing ultracold atom pair, such as in BEC, where both 
conditions 5 ~0 and a <C k c are well fulfilled, arbitrarily 
high atom-atom entanglement may be produced by a sin- 
gle photon scattering and detection, once the mismatch 
energy Em is well compensated [11] . 

For hotter atoms with temperature approaches or ex- 
ceeds the recoil temperature, or equivalent ly, cr > fc c , 
the scattered momentum wavepackets ^(q) will overlap 
each other and exhibit quantum interferences. From Eq. 
(fTT]h one sees that different orders of recoiled wavepack- 
ets take the same shape as the initial wavepacket G(q), 
together with an extra overall phase: (— l)' n l +1 Sgn(n). 
Explicitly, the first two orders of scattered wavepackets 
So(q) an d 5_i (q) take the same overall phase "—1", and 
therefore exhibit constructive interference; for higher- 
order wavepackets, however, S n (q) 1 take opposite phases 
"zbl" interchangeably, which will induce destructively in- 
terferences as in Fig. 3 (b). Due to this quantum inter- 
ference, the entanglement will be significantly decreased 
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FIG. 4: (Color online) The Schmidt number K for the en- 
tangled steady state of D(q, /c = k c ) is plotted against the 
reciprocal dipole parallelity parameter 1/5, with different ini- 
tial momentum widthes a as shown in the inset. 
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FIG. 5: (Color online) Time evolution of the produced en- 
tanglement. The dashed red line indicates E m /h — [llj ]. 
whereas the solid lines from top to bottom are plotted with 
£ m //i = 0.1xl0~ 3 (blue), £ m /ft = 0.5xl0~ 3 (green), E m /h = 
2xl0~ 3 (black), respectively. T = l, cr = 0.2/c c , 6 = 0.1. 



as in Fig. 4, therefore, one can hardly produce high en- 
tanglement for hot atoms with a ^> fc c , even when the 
conditions 5~0 and E m = [ll| are fulfilled. 

The time evolution of the produced atomic entan- 
glement can be obtained from the general solution in 
Eq. (|6|). As in Fig. 5, the steady entangled state is 
established with a time r which can be estimated as 



l/(£ 2 r). For h/E m ^>r [11], the maximal entangle- 
ment can be obtained by a single photon detection with 
the coupling time T satisfying r«T« h/E m . For the 
general case of Em ^ 0, the produced entanglement will 
first increase due to the pairwise scattering, and then 
decrease to a Bell-like state (K w 2) when T ^> h/Em, 
since the energy-conserved single-atom scattering pro- 
cess dominates in this time scale [9j . 

Conclusion. — From the first principle we have 
demonstrated that, the mechanism of pairwise photon- 



scattering in a couple of ultracold atoms may be used 
to produce superhigh atom-atom entanglement. When 
the atom pair is cooled well below the recoil temperature 
with parallel electric dipoles, it is possible to achieve an 
arbitrary amount of entanglement by a single photon de- 
tection once the scattering-induced energy mismatch E m 
is compensated [ll| . This scheme can be used to produce 
entangled atom resources for the test of EPR-nonlocality 
and for quantum information processing p], 0] . 

For experimental tests of this scheme, the spatially 
overlapping ultracold atom pairs can be produced by cou- 
pling weak atom laser beams from a BEC, and then be 
loaded into a cavity for detections of the scattered pho- 
ton [12|. Moreover, as in recent report [10[ of BEC su- 
perradiant 0, the generation of pairwise scattering by 
atom pairs is preferred over single atom scattering when 
the energy mismatch E m is compensated by the incident 
coupling light [10], therefore, it is most probable to an- 
alyze the quantum correlation in the pairwise scattering 
process with this proposed model. 
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